Download Room Acoustics Modelling using Gpu-Accelerated Finite Difference and Finite Volume Methods On a Face-Centered Cubic Grid In this paper, a room acoustics simulation using a finite difference approximation on a face-centered cubic (FCC) grid with finite volume impedance boundary conditions is presented. The finite difference scheme is accelerated on an Nvidia Tesla K20 graphics processing unit (GPU) using the CUDA programming language. A performance comparison is made between 27-point finite difference schemes on a cubic grid and the 13-point scheme on the FCC grid. It is shown that the FCC scheme runs faster on the Tesla K20 GPU and has less numerical dispersion than best 27-point schemes on the cubic grid. Implementation details are discussed.
Download Revisiting Implicit Finite Difference Schemes for Three-Dimensional Room Acoustics Simulations on GPU Implicit finite difference schemes for the 3-D wave equation using a 27-point stencil on the cubic grid are presented, for use in room acoustics modelling and artificial reverberation. The system of equations that arises from the implicit formulation is solved using the Jacobi iterative method. Numerical dispersion is analysed and computational efficiency is compared to second-order accurate 27-point explicit schemes. Timing results from GPU implementations demonstrate that the proposed algorithms scale over their explicit counterparts as expected: by a factor of M + 2, where M is a fixed number of Jacobi iterations (eight can be sufficient in single precision). Thus, the accuracy of the approximation can be improved over explicit counterparts with only a linear increase in computational costs, rather than the quartic (in operations) and cubic (in memory) increases incurred when oversampling the grid. These implicit schemes are advantageous in situations where less than 1% dispersion error is desired.
Download Non-Iterative Schemes for the Simulation of Nonlinear Audio Circuits In this work, a number of numerical schemes are presented in the
context of virtual-analog simulation. The schemes are linearlyimplicit in character, and hence directly solvable without iterative
methods. Schemes of increasing order of accuracy are constructed,
and convergence and stability conditions are proven formally. The
schemes are able to handle stiff problems very efficiently, because
of their fast update, and can be run at higher sample rates to reduce
aliasing. The cases of the diode clipper and ring modulator are
investigated in detail, including several numerical examples.